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Abstract 

The orbitally excited heavy quark baryons are studied in the Callan Klebanov bound 
state model with heavy spin symmetry. First, a compact description of the large iVc, 
infinite heavy quark mass bound state wavefunctions and the collective quantization is 
given. In order to study the kinematical corrections due to finite masses we motivate an 
approximate Schrodinger-like equation for the bound state. The effective potential in this 
equation is compared with the quadratic approximation (spherical harmonic oscillator) 
to it. This oscillator approximation is seen to be not very accurate. It is noted that the 
present experimental information cannot be even qualitatively understood with the usual 
light sector chiral Lagrangian containing only light pseudoscalar mesons. The addition of 
light vector mesons helps to overcome this problem. 
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1 Introduction 



The "bound state" picture [0, 0, in which a baryon containing a heavy quark is visuahzed 
as a bound state of a nucleon-as-Skyrme sohton with a heavy meson, is a very appeahng 
one. It has the interesting feature that experimental information from the mesonic sector 
of the theory (representing an approximation to QCD in the large limit) can be used 
to predict the properties of the heavy baryons. Originally it was applied Jl], |^ to studying 
the ordinary hyperons, but a straightforward extension was also made ^ to the c and b 
baryons. After the recognition of the importance of the Isgur-Wise heavy spin symmetry 
IQ, the study of the c and b baryons was pursued by seveal groups The relatively 

large number of papers suggests the richness and technical complexity of the approach. 



The present paper was stimulated by experimental evidence |]TT| for two orbitally ex- 
cited heavy c— baryons which may be interpreted as the heavy spin multiplet (A'^, A^*) with 
spin-parities (| , | ). Some properties were already studied in a bound state inspired 
framework based on a potential for the spatial wavefunction of the form 

V{r) = Vo + ^^r^ , (1.1) 

where Vq and k are constants. That treatment regarded Vq and k as arbitarary parameters 
to be fit. In the bound state approach they are, however, computable in terms of the 
Skyrme profile functions and the light meson - heavy meson coupling constants. While 
the latter especially are by no means precisely fixed it is easy to see that, in models 
in which the only light mesons present are the pseudoscalars, reasonable choices predict 
values for Vq and k which give an unbound or just barely bound ground state particle Ac. 
(The binding is actually sizeable: m{Ac) — m{D) — m{N) = —0.63 GeV). Hence there 
appears to be an important gap between the practical use of ( |1.1| ) and the actual bound 
state calculations. 

In order to investigate this problem we first formulate an approximate Schrodinger like 
equation which should hold for finite heavy meson mass M. The underlying equations 
that one gets, even for the ground state 0, are not of the simple Schrodinger form, but 
comprise three coupled differential equations. We therefore made a simple approximation 
which should be good for large M and which gives a Schrodinger- like equation with a 
potential function which may be compared with (|1 . 1| ) . It turns out that the quadratic 
approximation (|1.1| ) is roughly reasonable for the low-lying 6— baryon energies but is 
significantly worse for the low-lying c— baryon energies. The Skyrme potential gives more 
deeply bound baryons than (|1 . 1| ) ; however, the extra binding turns out to be nowhere 
enough to solve the problem. The non-relativistic form of the approximate Schrodinger 
equation enables us to easily make very important "two-body" corrections corresponding 
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to finite nucleoli mass (which is infinite in the large Nc starting point) by introducing the 
reduced mass. It is explained how this further reduces the ground state binding. The 
most straightforward way to obtain reasonable values of the binding energy seems to be 
to introduce light vector mesons |0, |[ in addition to light pseudoscalars. These provide 
a great deal of extra binding strength. The requirement of explaining the spectrum is 
noted to yield non-trivial constraints on the light vector meson sector of the theory. 

Another new aspect of this paper is the presentation of a somewhat streamlined ap- 
proach to the excited baryon wavef unctions. Following the approach of refs. 0, |^ for 
the ground state we show how excited state wavefunctions which are already diagonal 
(in contrast to those of |§ and may be written down almost by inspection. This is 
given in section 3 where the highly degenerate spectrum in the M —>■ oo, Nc —>■ oo limit is 
discussed, including the effects of light vector mesons. The physical states of the theory in 
this limit are recognized in section 4 after introducing collective variables to describe the 
Skyrme "tower" . This must be quantized as a boson so the low-lying states look like the 
quark model ones [1^ wherein a light diquark (belonging to the flavor SU{3) representa- 
tions 3 or 6) is rotating around a heavy quark with effective orbital angular momentum 
ieff- III turn, ieff equals the "light" part of the "grand spin" of the heavy meson in the 
background Skyrmion fleld. Finally, the approximate Schrodinger like equation for the 
physical kinematics situation is treated in section 5 and conclusions are drawn. 



2 Notation 

We will employ the same notations as in the earlier papers [0 and P]. For the reader's 
convenience, a very brief reminder is included here. 

The total effective chiral Lagrangian is the sum of a "light" part describing the three 
flavors u,d,s and a "heavy" part describing the heavy (0~, 1~) meson multiplet H and 
its interaction with the light sector: 

^eff flight ~l~ ^heavy • (^■-^) 

The relevant light flelds are the elements of the 3x3 matrix of pseudoscalars cj) and of 
the 3x3 matrix of vectors p^. Some objects which transform simply under the action of 
the chiral group are 

e = e'^/^- , f/ = , 

^M^(p) = d^p^ - d^Pf, - ig[p^, p^] , (2.2) 

where the pion decay constant F^^ ^ 0.132 GeV and the vector meson coupling constant 
g ~ 3.93. References on Cught may be traced from 0, |^. 



3 



The heavy muhiplet field combining the heavy pseudoscalar P' and the heavy vector 
Q'^, both moving with a fixed 4- velocity is given by 

H = (i^^)(z75P' + ^luQ'.) , S = 74ift74 ■ (2.3) 
In this convention H has the cannonical dimension one. For Cheavy we take: 

Cheavy ^ zV^TT[H{d^ - tagp^ - z(l - a)v^)H] + tdTT[H^^^,p^H] 



M 



tc — 
+—Ti[H^,j,F,,{p)H] , (2.4) 



where M is the mass of the heavy meson, rriy ~ 0.77 GeV is the light vector meson mass 
and 

v,,p,= '-i^d,^^±^^d,0 ■ (2.5) 

d, c and (ag) are respectively dimensionless coupling constants for the if— light pseu- 
doscalar, if— light vector magnetic type and if- light vector gauge type interactions. It 
seems fair to say that these coupling constants are not yet precisely fixed from experiment. 
For definiteness we shall use the values 

rf = 0.53, c=1.6, a = l. (2.6) 



The values for d and c are based on single pole fits |14, 15] to the experimental data [16| on 
the D ^ K and D K* semi-leptonic transitions. The ratio c/d is consistent with that 
obtained [T^ using a suitable notion of light vector meson dominance for the D* D'j 
and D* Dtt branching ratios. It should be remarked, however, that some authors have 



suggested [jTH| a smaller value of d, around 0.3. The value a = 1 is based on light vector 



meson dominance which seems reasonable. Earlier we used a negative a based on 
fitting the heavy baryon properties without using the excitation energy constraint. The 
present paper contains a reconsideration of that fit. 



3 Excited Baryon States at the Classical Level 

In this section we first write the classical Skyrme solutions for the light part of the action 
(which is taken to include the first three fiavors). Then we find the wavefunctions for which 
this "baryon as soliton" is bound to the heavy meson H. All the orbitally excited states, 
rather than just the ground state, will be included. We shall work in the heavy quark 
symmetry limit so that the radial wavefunctions reduce to their delta-function limits. It 



4 



turns out that the bound excited wavefunctions are remarkably simple generalizations of 
the ground state one given in and (see section 3 of 0, for example). 
The usual hedgehog ansatz for the light pseudoscalars is 

j^exp[.x^r^] Oj _ ^3_^^ 

When we include the light vectors, we have similarly the classical solutions 

ft..= (**"-0^"°''-' ?) . (3.2) 

with 

Pic = /l~ (^ikaXkG{r) , Poc = 0^ iLiic = 0, uJoc = (^{r) . (3.3) 
V2gr 

The appropriate boundary conditions are 

F(0) = -7r, G{0) = 2, u;'(0) = 0, 

F{oo) = G{oo) = uj{oo) = . (3.4) 

Now, following the Callan-Klebanov approach, we want to find the wavefunctions of a 
Schro dinger- like equation for H in the classical background field above. Since the above 
ansatzae mix isospin with orbital angular momentum, it is very convenient to make a 
partial wave analysis in terms of the grand spin G, 

G = I + L + S, (3.5) 

where L is the orbital angular momentum of the heavy meson field, S is its spin and 
I its isospin. Due to the heavy spin symmetry, the portion of S from the heavy quark 
decouples from the problem. Thus with the decomposition 

S = S' + S" , (3.6) 

where S" is due to the heavy quark while S' is the remainder, the truely relevant object 
is the "light grand spin" : 

g = I + L + S'. (3.7) 

We pointed out earlier [0 that the ground state wavefunction is characterized hj g = 0. 
Now we want to study the excited states. Remember that in the heavy meson rest frame 
the 4x4 matrix H has non-vanishing elements only in the lower left 2x2 sub-block: 
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Here the index / represents the spin of the hght degrees of freedom within the heavy 
meson, the index h stands for the heavy quark spin and the index h represents the isospin 
of the hght degrees of freedom within the heavy meson. FoUowing |^ and |p we write the 
general 7^ wavef unction as: 

ft) = I ® ■ ■ " = f (3.9) 

wherein the radial function u{r) satisfies r'^\u{r)\^ 5{r) and Xh is the heavy quark 
spinor. We will see how the presence of the x ■ r factor simplifies the calculations. The 
remaining factor ip^i is a kind of generalized spherical harmonic with the appropriate 
covariance properties. There is the standard three-fold ambiguity of which two vectors in 
(|3.7| ) should be coupled together first. Here, coupling together 



K = I + S' (3.10) 

leads to a nice further simplification. Then we write, for a ipdi with orbital angular 
momentum r and "K-spin" k in ( p.lO| ) coupled to light grand spin g: 

gs; r,k) = J2 Cl^ks^g-Jr^viidAk, k^) , (3.11) 

where C stands for the Clebsch-Gordon coefficients, Y stands for the usual normalized 
spherical harmonics and ^di{ki k^) are the isospin-light spin wavef unctions: 

1 

0) = ^{sy, + s\tl), -1) = s\ii , (3.12) 



e,KO, 0) = ^{s]i\, - s\i\), 1) = s]i\ , 



s and i being the two component S' and I spinors. Notice that multiplying tjjdi by (x- T)arf, 
as in (|3.9| ), does not change its (?— spin. 

So far, we have just written down possible wave functions allowed by (7— spin covari- 
ance; the question of which ones or linear combinations actually correspond to bound 
states has not yet been addressed. To investigate this question we consider the matrix 
elements of the "potential" V obtained from Cheavy in ( P-4| ). (See (3.2) and (3.3) of [0 for 
further details). The potential conserves (7— spin and parity so we need only its matrix 
elements between wavefunctions H and H' of the same g and parity: 

-MdF'{0) 



dF'{0) 



J d^X {HXi' (x ■ T)a'b (Tin ■ Tbc (x ■ T)ca {H)t^ + ■■■ 

dn {^%, an ■ Tbc (^)az (X'^X) + ■ ■ ■ , (3.13) 
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wherein dO, denotes the angular integration and the second hne follows from the first 
on the substitution of (|3.9|) and the observation that the definition H = 74-^^74 for the 
4-component matrix H results in an extra minus sign for the two component if^;. The 
three dots in ( p.l3| ) stand for contributions arising from the couplings of the light vectors 
to the heavy meson; these do not change our conclusion and will be discussed later. 

Notice that all the x ■ r factors have disappeared in the second line of ( ^.131 ). This 
means that the effective potential operator for the wavefunction (|3.11|) is simply: 



irfF'(O) o- • T = idF'(O) [2K2 - 3] , (3.14) 

where we have used the fact that cr ■ r is acting on wavefunctions of definite K. For A; = 
we have the energy eigenvalue — |(iF'(0) ~ —0.63 GeV, which indicates that the k = 
states are the bound ones. The (in general) three k = 1 states (corresponding to r in 
( |3.11| ) taking on the values g — 1, g, and g + 1 ) have the energy eigenvalue +^gF'{0) and 



are unbound in this model. The wavefunctions ^pdl in (|3.11D as well as the if's in 



are already diagonal - a circumstance following from the choice of k as the intermediate 
label for constructing states of good g. It is not actually necessary to carry out the 
multiplication by x ■ r in (|3.11| ); however, this is done in Appendix A. We shall refer to 



the "fc— value" of a wavefunction as that of the factor ipdi in ( p.ll|) . As seen in Appendix 
A, this gets modified on multiplication by x ■ r. 

The bound state wavefunctions may be written in a very simple form. Since they have 
k = 0, ( |3.11| ) simplifies to Ygg^ ^^li^y^)- Furthermore ^d«(0)0) = "TS^"^'' S^^i^S finally for 



dp) 



u{r) 



X ■ T)ad ^dl yggs Xh , C = 1, 2 



, a = 3 . 



The g = ground state wave function is seen, using loo = l/v47r, to coincide with (3.7) of 
@. In fact, the bound orbitally excited wavefunctions are simply obtained by multiplying 
the ground state one by \/A7iYgg^. 

The parities of the bound state wavefunctions are given by the formula: 

parity = (-1)^ . (3.16) 

This follows most directly from the fact that a negative parity meson is being bound in a 
linear combination of states with orbital angular momenta g — 1 and g + I (see the first 
line of (Al)). Of the three (in general) unbound states with given g, one has parity (— 1)^ 
and the other two have parity —(—1)^. 

The light grand spin quantum number g has yet further physical significance for the 
bound state wavefunctions. Note that in the extreme limit in which we are presently 
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working all the bound states are degenerate in energy. This degeneracy will be broken if 
one allows for finite masses. We then expect a centrifugal contribution to the energy of 
the form 

^4//(4// + l), (3.17) 

which would raise the energies for orbital excitations. Callan and Klebanov pointed out 
in their original paper [||] that £e// differs from the orbital angular momentum due to the 
isospin-angular momentum mixing in the Skyrme ansatz. It turns out that, in fact, 

ieff = 9. (3.18) 

To see this, note that the centrifugal energy operator is approximately at small r given 

by 

1 



+ 2 + 41 ■ L) 



(3.19) 



2Mr2 

as is mentioned in [|l| and as also emerges in the approximation for finite M which perserves 
the heavy quark spin multiplets (to be discussed in Section 5 of this paper). With A = L+ 
I, (|3.19|) may be rewritten as (2 — + i)/ (2Mr'^). Now the bound state wavefunction 
is a linear combination of i = g — 1 and i = g + 1 pieces. Considering a basis in which A 
is diagonal, we note that the i = g — 1 piece can only couple to light grand spin g for the 
choice X = g — ^. Then the centrifugal energy becomes 



1 



'2i9-l)i9 + l)-i9-'^)9 + l 



2Mr 



^9(9 + 1) 



(3.20) 



g + \, which leads to the same result. 



2Mr2 

Similarly the i = g + 1 piece requires A 

To end this section we give the form of the potential operator when the effects of light 
vector mesons are included. Then ( |3.14| ) should be replaced by: 



cr ■ T 



-dF'iO) 



rriyg 



+ 1 



^9 



(3,2i; 



In this formula d, c and a are the heavy meson— light meson coupling constants defined in 
(|2.4| ), while F{r), G{r) and uj{r) are the pseudoscalar, p meson and u meson soliton profile 



functions defined in (|3.1| ) — ( p^ . Notice that the p meson piece has the same cr ■ r factor 
as the pseudoscalar piece while the uj meson piece has simply a 1 in the light spin— isospin 



spaces. This shows that the wavefunctions of (|3.9|) and ( 3.11 ) are still diagonal. The 
eigenvalues for the = states and for the k = 1 states are now read off to be: 



V{k = 0) 
V{k = 1) 



(3.22) 
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The quantitites F'{0), G"{0) and ci;(0) are obtained by solving the coupled differential 
equations which arise by minimizing the static energy of Cught (describing the ordinary 
baryons). This yields for a typical best fit to baryon and meson masses: 

F'(0) = 0.795 GeV, G"{0) = -0.390 GeV, uj{0) = -0.094 GeV . (3.23) 



It should be remarked ||2^ that the sign of a;(0) is linked to the sign of the puj(f) 
coupling constant; the sign has been chosen EDI to yield a best fit to light baryon 



electromagnetic form factors. To be on the conservative side one might want to consider 
the possibility of reversing this sign. Now, substituting (|2.6|) and ( 3.23 ) into ( 3.22| ) yields 



V{k = 0) = -0.63 - 0.62 + 0.26 (-0.26) 
V{k = l) = +0.31 + 0.21 + 0.26 (-0.26) 



-0.99 (-1.51) GeV, 
-0.78 (0.26) GeV, 



(3.24) 



wherein the order of the terms in ( p.22|) has been retained. The numbers in parentheses 



correspond to reversing the sign of uj{r) in |3.23| . We may observe that the p meson 
contribution (G"'(0) terms) strengthens the attraction in the bound channel and also 
increases the repulsion in the unbound channel. The u term is not dominant but may 
have significant effects. The binding in the present approximation seems somewhat too 
large. However, finite M effects will provide a substantial reduction. 



4 Physical States of the Model 

The states of definite angular momentum and isospin emerge, in the soliton approach, 
after collective quantization. The collective angle- type variable A(t) |^l], 22] is introduced 
on the light fields as 



e(x, t) = A{t)U^)A\t), p^(x, t) = A{t)p^,{^)A\t) , 



(4.1) 



where Cc(x) and p^d^) are given in (|3.1|) - (|3.3|) and generalized angular velocities flk are 
defined by 



(4.2) 



k=l 



the Afc being SU{3) Gell-Mann matrices. All the heavy meson states, rather than just the 
ground state, may be included by introducing a Fock representation for the field H as 



^(x,t)=^^„(x)e^^"*at 



(4.3) 
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wherein the stationary states n = {g,g3;r,k,s'^} are given in ( |3.9| ) and (the un- 

bound states are also included) while aj^ is the creation operator for the state n with 
energy eigenvalue En- The "light part" of H is subjected to a collective rotation with the 
replacement 

H^^ii^,t)^Aat{t)HU^,t). (4.4) 



Now the collective Lagrangian, Lcoii, is obtained by substituting (|4.1|) and ( [4 .41) into 
/ cP X {Ciig fit + 'C heavy) and carrying out the spatial integration. The final physical states are 
recognized after the quantization of Lcou- There are two ways in which this Lcou differs 
from that of the usual SU{3) Skyrme model |Q. The first is that instead of a term 
— ^ ^8 there is now (in the one heavy quark subspace) a term — ^ i^g- As discussed in 
section 4 of , the physical significance of this fact is that there now exists a constraint on 
the quantum states which requires the "Skyrmion rotator" to transform as an irreducible 
representation {^} of SU{3) which contains a state with hypercharge Y = 2/3. This 
particular state must necessarily have integer isospin. Furthermore its spin, denoted J^, 
must equal its isospin. The two lowest Skyrmion multiplets (higher ones are probably 
model artifacts) are 

(/i = 3, Js = 0) 

(/i = 6, J, = 1). (4.5) 

Note that the Skyrmion rotator behaves as a boson; at the collective level it evidently 
describes a diquark state. This analysis is forced upon us when we consider the case of 
three light flavors with the attendant Wess-Zumino term. It agrees with the original 
picture of "spin-isospin" transmutation; namely, at the collective level the heavy meson 
field 

(i) loses its flavor quantum numbers 

(ii) acquires a spin equal to G (4.6) 



The full baryon state is a product of the bosonic diquark rotator in (|4.5| ) and the fermionic 
bound state wavefunction corresponding to ([4.6| ). 

The second way in which Lcoii differs from that of the SU{3) Skyrme model is the 
presence of the following term linear in Q: 



r,(fi) = (^-a)fi, -(l-a)x-f7% , (4.7) 
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where H is taken from ( [4 .31) and rj(f2) is evaluated near the origin of x-space. From the 
analysis of one might expect a term like ([4 .71) to lead to (hyperfine) splitting of the 
baryons belonging to a given heavy spin multiplet. Fortunately, this does not happen 
as one may see by inspection using the bound (negative energy) wavefunctions of ( p.l5| ). 
Then ( [4. 7] ) contains a factor: 

ed'i (x ■ T)d'a ir^)ab (x ■ t)m = Tr(x ■rr^x-r) = 0. (4.8) 

Thus, in the heavy mass M oo limit, the bound eigenstates of the collective Hamil- 
tonian are simply products of the wavefunctions ( |3.15 ) with the di quark Skyrme rotator 



wavefunctions. Taking account of ( [4.6| ii), the total angular momentum J is given by 

J = J, + G = J, + g + S", (4.9) 

where, for convenience, the heavy quark spin S" has been made explicit in the last step. 
The Skyrme rotator wavefunctions are p4|: 



^rot{fi,YIh,JsM) = i-iy-'-'^dtrnfi D'^'*^^,, j^ _^^iA), (4.10) 

where D^^\A) is the representation matrix of SU{3). Defining the total light spin j = 
J — S", we write the overall wavefunctions for the bound states as: 

**(/x;F,J,J3; J.,^,j;j3,4') = E Cj^l'^^j^Kotif^^YIh, JsM) H^.ig, g,, s'^ , (4.11) 

where Hj^f^ is given in ( |3.15| ). Finally j and S" may be added using C/^^ss-.Ja to yield the 
heavy spin baryon multiplets having J = j ± |. 

All the states in ( [4.11[ ) have, up to relatively small 0{1/Nc) corrections, the degenerate 



energy eigenvalues V{k = 0) in (|3.22|) . The 0{1/Nc) corrections associated with the 
present collective quantization splits the 3 and 6 SU{3) representation states according 
toi: 

M{6) - M{3) = -[m{A) -m{N)] , (4.12) 
3 

where A and stand for the A(1230) and nucleon masses. Now let us ennumerate the 
physical states. We expect, as mentioned after (|3.17| ), that finite M will split the huge 



degeneracy so that g = corresponds to the ground states, g = 1 to the first excited 
states and so on. For g = the Clebsch Gordon addition in ( |4.11| ) becomes trivial and 
the discussion reduces to that given in section 4 of All the g = states have (see 
( |3.16| )) positive parity. The SU{3) representation 3 with j = Js = has the content 
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{Aq, Sq(3)}. It has spin J = ^ obtained by adding the heavy spin s" = | to j = 0. 
The SU{3) representation 6 has j = Js = I. Adding s" to this yields a degenerate heavy 
spin multiplet with both J = \ and ^ = | members. The flavor content is denoted 
{Sq, Hq(6), Qq} and {S^, 2^(6), Q*q}. 

It is sufficient to consider the ffist orbitally excited states with g = 1 in order to 
see the general pattern. All the g = I states have, according to ( p.l6| ), negative parity. 
The SU{3) 3 diquark rotator state with = is coupled in ( |4.11| ) to total light spin 
j = g = 1. Combining this, in turn, with the heavy quark spin S" yields a heavy spin 
multiplet with both J = ^ and ^ = | members. These are denoted as {Ag, Sq(3)} and 
{A'^, S'^(3)}. The SU{3) 6 diquark rotator state with = 1 is coupled in ( ^TlD to total 
light spin j = 0,1 or 2. Combining each of these three with S" yields three heavy spin 
multiplets with total spin contents (J = |), (J = |, J = |) or (J = |, J = |). For 
each of these there are 6 flavor states analogous to {Sq, Sq{6), ^q}. Notice that the 



zero isospin states denoted A'^ and A'^* are possible candidates for recently discovered fTT 
resonances. 

SU{3) splittings in the M — ^ cxd limit have been discussed in section 5 of for the 
g = states; a similar analysis can be given for the orbitally excited states. 

The positive energy (unbound) states, corresponding to A; = 1 in ( |3.11| ), can be ennu- 
merated in a similar way using the appropriate Hf^{g, g3;r, k, s'^) in ( l.llj ). While these 



states would presumably not be bound in the true theory, they can be expected to play 
a role as virtual intermediate states for further application of the formalism. 

5 Kinematic Corrections 

The preceeding results hold in the large Nc limit (where the nucleon mass is formally 
infinite) and in the large M limit (where the heavy meson mass is formally infinite). In 
order to compare with experiment it is clearly important to get an idea of the corrections to 
be expected in the real world. A complete discussion of these effects would be enormously 
complicated and beyond the scope of the present investigation. Thus we will content 
ourselves with a somewhat schematic model which has the advantage of simplicity and 
which is familiar enough to stimulate our intuition. The most straightforward way to 
proceed is to start with the simplest "ordinary field" Lagrangian which reduces to the 
heavy field Lagrangian ( |2.4|) in the M ^ oo limit. Such an ordinary field Lagrangian, 
constructed with a heavy pseudoscalar SU (3) triplet P and a heavy vector SU (3) triplet 
Qn was given in (3.25) of |T^ in connection with our original discussion of ( |2.4| ). For 



additional simplicity we shall at first neglect the heavy meson interactions with the light 
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vectors, although we expect from (|3.22|) and (p.24|) that they are actually important. 



Then the heavy Lagrangian becomes: 
Hheavy = 'V^PV^P - M^PP - ^ {V^Q, - V,Q^) (v^Q, - V,Q^) - M*^Q^Q^ 
+2iMd (Pp^,Q^ - Qf,Pf,P) - id'ef^apt, ippQaPt^Qp - QaPpT^pQp) , (5.1) 

where "D^P = d^P — iv^P, etc. Note that and p^ are defined in ( p.5| ). In writing ( |5.1|) 
we have allowed for different P and masses M and M* as well as different coupling 
constants d and d' . But for simplicity we shall further restrict M = M* and d = d'; thus 
our model will treat only those corrections to the heavy quark symmetry due to finite M. 
The same model ( p.l| ) has been investigated by Oh et. al. |^ for the ground state heavy 
baryons and with choices oi d ^ d' and M ^ M' made to fit the heavy baryon masses with 
experiment. Our treatment will differ in a number of ways. First we shall approximate the 
three coupled differential equations which result by a simple Schro dinger- like equation. 
We shall generalize the approximate equation to include excited states, and compare with 



the more standard picture |12| of a quadratic potential. We will also give some discussion 
of corrections due to the finite nucleon mass. 

To begin, we look for the stationary ground state solutions of the equations of motion 
(see Appendix B) which result from ( ^.1|) . We know what this should look like in the 
heavy limit. Remember that the field P above is related to the field P' in (|2.3| ) by 

P = e^'^^-^P' , (5.2) 

and similarly for Q^. The ground state {g = 0) wavefunction for H in (|3.15|) translates to 

Qf = — Tr(7,^) = -^4i= [x, 5m ~ ^(x x r,,),] p, , (5.3) 
2 2 VSvrM 



where b is the 2-valued isospin index and pd = edeXe- Taking ( p.2| ) into account and 
allowing different radial dependences for the different terms in ( |5.3| ) suggests the ansatzae 
n for the ordinary (unprimed) field stationary solutions: 



pb _ VV) _\ ^ Auit 

r — —=[:S.- TjhdPdG 

v47r 



Pde'^\ (5.4) 



'Air 

where u is the relativistic energy. At the ordinary field level we should interpret pd as the 
isospace spinor. With this parametrization the M — > oo limit is expressed as 

0=-V^i = -^. (5.5) 
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Substituting 
equation: 

<p +- 



into the P equation of motion given in (B.l) yields the differential 



+ 



y/2MdsmF 



1 r 1. 



= 



cosF)(3 + cosF)} 



- MdF'ipi 



(5.6) 



where a prime denotes differentiation with respect to the radial coordinate r. The two 
other similar coupled equations which result from (B.2) are given in (B.4) and (B.5). First 
let us check the consistency of the M — * oo limit. Dropping all terms not proportional to 
M and substituting ( ^Tsf ) we easily see that each of (|5.6| ), (B.4) and (B.5) reduces to 

2sinF 



(a;' - M') + Md{F' - 







(5.7) 



This equation has a solution where 0(r) is sharply peaked around the origin and where 



E, 



2 



F' - 



2 sin F 



Here we have defined 



E. 



M2 



e// 



2M 



r=0 



LO 



--dF'iO) . 



M. 



(5. 



(5.9) 



(|5.8| ) is seen to agree with the binding energy V{k = 0) in (|3.22|) , when the contributions 
from the light vectors are neglected. 

Now let us approximate ( p.6|) to achieve an easy form for convenient further study. 
It seems very natural to isolate the effects of finite M by retaining it in ( p.6| ) while 
approximating ipi and tp2 by their M oo forms from ( |5.5| ). Then we end up with a 
completely standard looking Schrodinger equation for the radial wavefunction. We set 
0(r) = w{r)/r as usual, divide through by 2M and finally add the "centrifugal" term 
( |3.17| ) so as to generalize the equation to include orbitally excited partial waves. Then 
( |5.6| ) becomes 



2M 



w"{r) + 



4//(4// + l) 
2Mr2 



eff\ 



w{r) = Eeffw{r) 



(5.10) 



wherein, 



eff\ 



-d 



F' - 



2sinF 



Mr2 



1- -(1-cosF) (3 + cosF) 



(5.11) 



Of course, w{r) has an implicit ^eff label. Even though ( ^.10|) looks exactly like a non- 
relativistic Schrodinger equation it actually (noting ( ^.9|) ) contains the energy in the 
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characterestic relativistic manner, {uJ^ — M^). It describes a meson of mass M in the 
potential field V^//(r) due to an infinitely heavy (large limit) nucleon. Notice that the 
piece inVeff which has an overall ^ factor behaves as for small r; this term is in 
fact very small. There is no exact analytic form for the Skyrme profile F{r). To make 
the equation self-contained we adopt the Atiyah-Manton approximation p5|: 



F{r) — TT 



(A2 + r2)i/2 



(5.12) 



where the parameter choice = 15.61584 GeV ^ corresponds to F'ifS) = 0.795 GeV from 
( |3.23| ). Using (|5.12|) , the effective potential V^//(r) is graphed in Fig. 1. Also shown is 



the quadratic approximation which leads to a spherical harmonic oscillator potential: 

VsHO = -ldF'{0) + ^Kr\ (5.13) 



where k is given in (4.4) of [|7|. With the neglect of the light vectors and the choice ( p.23| ), 
we have K = 0.1562 GeVl It is seen that the quadratic approximation is not a very 
accurate representation away from small r. This region is especially relevant for orbitally 
excited states. On the other hand we note that, as expected, the Ve//(r) which results 
from the Skyrme approach is not confining. Thus it is probably not trustworthy for higher 
orbitally excited states. 

To avoid confusion we remind the reader that, according to the discussion of section 
4, there are many physical states which correspond to a given ieff- They all have parity 
(^—lyeff^ 'j]2g heavy baryons which belong to the SU{3) 3 representations comprise a 
heavy spin multiplet with spin content (4// ~ |)) (4// + |) [for 4// = this collapses 
to just |]. The heavy baryons which belong to the SU{3) 6 representation form three 
heavy spin multiplets with spin contents [{ieff — f), (4// — I)], [(4// ~ |)? (4// + |)] 
and [{ieff + {ieff + §)] [for ieff = this collapses to just (|, |)]. 

Now let us consider the numerical solutions of ( p.lOp - ( |5.12| ). For definiteness we take 



the weighted average of the pseudoscalar and vector meson masses to obtain M; this 
yields M = 1.94 GeV for the charmed meson mass and M = 5.314 GeV for the 6— meson 
mass. The wavefunctions w{r) are taken to behave as A^'^ff'^^^ at the origin. Displayed in 
Table 1 are the values of Eeff corresponding to the most deeply bound state in each of the 
ieff = 0, 1, 2 channels for both the c— baryons and 6— baryons. For comparison we also 
show the results obtained from the spherical harmonic oscillator (SHO) approximation 
(|5.13| ). In that case we have a well known analytic formula for the energy levels: 



^{N) 3 . rK~ /3 



Si;; = --dF'(0) + ^- - + . (5.14) 
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4// 


Ee// from (Km 


SHO approximation 




c — n fi T\TC\X\ t; 

o UfXl yvjlio 




n 


— n 977 


—0 904 


1 


-0.113 


+0.079 


2 


-0.012 


+0.363 




6— baryons 







-0.403 


-0.373 


1 


-0.277 


-0.205 


2 


-0.169 


-0.031 



Table 1: Eeff in GeV from (5.10) and from the spherical harmonic oscillator (SHO) 
approximation. 



where the level corresponds here to states with if;ff = N,N — 2,---. It is interesting 
to note that the exact numerical ground state {ieff = 0) energy eigenvalues are more 
negative than the approximate SHO ones. This can be physically understood from Fig. 
1, since the SHO potential is narrower and should thus have more zero point quantum 
fluctuation energy. As we expect, the SHO approximation is better for the 6— baryon 
case than for the lighter c— baryon case where the more spread out wavefunctions sample 
more of the differing large r regions. Similarly, the SHO approximation is better for the 
ground state than for the excited states. In fact the c— baryon i^ff = 1 state is unbound 
in the SHO approximation. The model also predicts radially excited wavefunctions. For 
the c— baryon another £e// = level is found at —0.051 GeV. In the SHO approximation 
this is expected to be degenerate with the unbound ieff = 2 level at +0.363 GeV. As 
M — >• oo the SHO approximation results and the Skyrme potential results get closer to 
each other. For example, already at M = 15 GeV the difference between the predicted 
values of the ground state Eeff differ only by 1.5%. However, the collapse to the M ^ oo 
value of —0.63 GeV is rather slow; even at M = 100 GeV, there is a 10% difference. Note 
that the finite M effects are large both for the c and b baryons. 

A crucial question, of course, is how well these results agree with experiment. For this 
purpose it will be seen to be sufficient to neglect the relatively small 1/Nc corrections (of 
order 0.1 GeV). Then the data at present yields just three relevant numbers. First there 
are the "binding energies" 

B.E.c = m{Ac) -m{N) -m{D) = -0.63 GeV , 

B.E.b = m(Afe) -m(A^) -m(5) = -0.78 GeV. (5.15) 

In addition, if the recently discovered [|1^ zero isospin heavy baryons are identified with 
the A^ and A^* mentioned in section 4, we have the excitation energy 



E.E.e = m{A'^) - m(Ae) = 0.31 GeV . 



(5.16) 
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It seems natural to identify the binding energy with E^ff. Then the comparison with 
experiment is presented in Table 2. In column (A) the M — > oo results are given. The 
agreement with the binding energies is reasonably good but the excitation energies are 
predicted to be zero. The results of taking finite M corrections into account via the 
numerical solution of ( p.l(J| )-( p.l2| ) are shown in column (B). We note that the general 
trend of the experimental data is reproduced in the sense that each prediction of the 
model is about half of the experimental value. The experimentally greater binding for the 
b baryons compared to the c baryons is reproduced. This picture is very suggestive since 
we have been, for simplicity, working in a model which does not include the light vector 
mesons. Equations ( 3.22| ) and ( 3.24|) show that one may expect the binding to be greatly 
strengthened by the addition of the light vectors. To test this roughly one may use the 
SHO approximation including light vectors. Taking the inputs from (2^) and (|3.23| ) we 
have, in addition to V{k = 0) = -0.99 (-1.51) GeV from ( ^1241) , k = 0.295 (0.375) GeV^ 
from (4.4) of 0, resulting in 



^2 



VsHO (with light vectors) = -0.99 (-1.51) + 0.295 (0.375)— , (5.17) 

where the numbers in parentheses correspond to the reversed sign for ti;(r). (Of course it 
is equivalent, from the present standpoint, to keep the sign of uj{r) and reverse the sign 
of a.) 

The predictions from this model are shown in column (C) and are seen to be signifi- 
cantly closer to experiment. (The values in parentheses show the sensitivity to changing 
the sign of ci;(r).) It thus appears that the Callan Klebanov approach including light 
vector mesons and corrections can roughly explain the general features of the presently 
observed heavy baryon mass spectrum. This picture is however achieved in the ^ oo 
limit in which the nucleon mass is formally infinite. From a naive kinematical point of 
view this seems peculiar, although we may perhaps argue that the Nc oo limit is often 
more accurate than it has a right to be. 

It is clearly of interest to explore the "two body" corrections corresponding to taking 
the nucleon mass to be its finite experimental value. The most straightforward way to 
proceed is to note that ( [5.10D may be regarded as a non-relativistic Schrodinger equation 
with Eeff the binding energy. Then we replace M in ( |5.m| ) with the reduced mass 
/i = {Mm{N))/{M + m(A^)) according to the usual prescription. If this is applied to 
finding the ground state £e// = energy for the heavy charmed baryon (where fi = 0.633 
GeV) we find E^jf = -0.084 GeV rather than the value -0.277 GeV listed in Table 1. This 
drastic reduction in binding strength can be easily understood in the SHO approximation. 
There the classical binding energy of —0.63 GeV is pushed up by the zero point energy 
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Expt. 


(A) 


(B) 


(C) 


(D) 


B.E.c 


-0.63 


-0.63 


-0.277 


-0.41 (-0.85) 


+0.31 (-0.36) 


B.E., 


-0.78 


-0.63 


-0.403 


-0.64 (-1.11) 


-0.08 (-0.48) 


E.E.c 


0.31 





0.164 


0.39 (0.44) 


0.68 (0.77) 


E.E.b 


? 





0.126 


0.24 (0.27) 


0.61 (0.69) 



Table 2: Comparison with experiment in the large Nc limit. All quantities are in GeV. 
Column (A) gives the prediction in the M oo limit, column (B) corresponds to the 
solutions of (5.10)-(5.12). Column (C) corresponds to the SHO approximation when light 
vectors are included. Column (D) further includes Nc subleading two-body corrections 
due to finite nucleon mass to the approximation (C). 



However the replacement M fi greatly increases the zero point energy. 



To get a somewhat more realistic impression we use the SHO approximation with 
inclusion of the light vectors (|5.171 ). Then we obtain for the binding and excitation energies 
the results in column (D) of Table 2. It is clear that choosing the solution in parentheses, 
corresponding to the choice of positive u;(0) in (|3.23|) , yields a fairly reasonable picture 
for the ground state energies. The binding energies are somewhat too small in magnitude 
but we have seen already for the model without light vectors (Table 1) that there is 
more binding for the c— baryons with the Skyrme profile than is indicated by the SHO 
approximation. Thus we expect some improvements when we consider, in the future, 
the analog of ( |5.10| ) with the inclusion of light vectors in Ve//(r). Similarly, we expect 
considerable reductions in the predicted excitation energies. More precise determinations 
of the coupling constants d, c and a is evidently also important. Furthermore, relativistic 
two body effects may play a role. 

It does not seem that an adequate description of the spectrum can be obtained in 
a model with just the light pseudoscalars included and when finite m{N) is taken into 
account - the lowest lying charmed baryon would be barely bound. 

In this paper we have presented a streamlined formalism for investigating, in the Callan 
Klebanov bound state picture, the orbitally excited baryons containing a heavy quark. 
The effects of light vector mesons and strange quarks were included. Wavefunctions which 
diagonalize the bound state Hamiltonian in the M oo, Nc oo limit were obtained in 
a transparent way. These correspond to binding energies which are somewhat too large. 
Furthermore there is no splitting between any excited states and the ground state. For the 
purpose of investigating the more realistic case we developed an approximate Schrodinger 
like equation to describe the finite M situation, but where there is no breaking of the 
heavy spin multiplets. The significant difference between the potential function of this 
equation obtained from the Skyrme profiles of the light pseudoscalars and the spherical 
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harmonic oscillator approximation to it was pointed out. This was more important for 
the c— baryons and for the excited states. The SHO approximation may be considered 
reasonable only in a rough qualitative sense. The general effect of going to finite M, 
but still infinite m{N), was to reduce the strength of the binding. It was seen that the 
model including only light pseudoscalars was unable to give sufficient binding. With 
light vector mesons, however, the binding was substantially increased. At this level the 
excitation energy was predicted to be somewhat too low. Next, the effect of finite m{N) 
was taken into account in a non-relativistic approximation. This further decreased the 
binding strength and favored one variant of the model with light vectors. 

In subsequent work we plan to study the effects of using the Skyrme rather than the 
SHO potential (with light vectors) on the spectrum as well as on the wavefunctions (and 
Isgur-Wise functions). We would also like to further study heavy spin and SU (3) breaking 
in the present framework. 

Appendix A 

Here we display the effect of multiplying the wavefunctions ipdi{g,g3',r,k) (defined in 
( |3.11| )) by (x-T)arf, as required for the overall wavefunction given in ( p.9| ). For simplicity we 
shall restrict ourselves to the case gs = g and use the shorthand notation i^diig, gs', k) = 
$(r, k). Then 



(x-t)$(^7,0) 
{±-T)^g,l) 
(x-t) $((7-1,1) 
(x-t)$((7 + 1,1) 



V2g + T 

-1 

V^gTT 
1 

V^gTT 



y^<!?{g-l,l)-^g + l<!?ig + l,l) 
7n$((?-l,l) + v^<l>((7 + l,l) 
gTl^g,l)-,/^Hg,0) 



-^^g,l) + Jg + l^g,0) 



(A.l) 



The wavefunctions in the first line, with A; = 0, are the bound ones. The others, with 
k = 1, are unbound. Note that the wavefunctions on the first two lines both have 
parity = (—1)^ and are manifestly orthogonal to each other. Those on third and fourth 
lines have parity = —(—1)^ and are also manifestly orthogonal. In the special case when 
g = there are just two, rather than four, independent states - $(1, 1) with positive 
parity which is bound, and $(0,0) with negative parity which is unbound. 

Appendix B 

The equations of motion which follow cannonically from ( p.l| ) are: 

-V^V^P + M^P-2iMdp^Q^ = 0, (B.l) 
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+2id' epy^^ Pf,T>^Qi3 = . 



iB.2] 



We choose to deal with P rather than P since P transforms hke a hght quark (rather 
than antiquark) and the interesting heavy meson dynamics is associated with its hght 
constituents. (|B.2|) does not involve time derivatives of the component (54- Hence is 



non-dynamical and may be eliminated in terms of the other fields. Correct to subleading 
1/M order, 

1 



Q4 



M 



*2 



[B.3) 



Making the ansatzae ( ^.4| ) for the ground state wavefunction and substituting into ( |B.1| )- 
rsD yields, in addition to (|5^) the two additional equations: 



-^'1 - + 
r 



M 



*2 



+ +sin^^l 



+V2 



2 sin^ f 



+ 



F' sin F ud' sin F 



2r 



+ 



^1 

iIj2 + MdF' i 



(B.4) 



r 



M*2 - cu^ + 4(1 + 2 sin^ - - 2 sin^ -) - cud'F' 



^2 



+V2 



2sin2| F'sinF ood' sinF 
^ + — + 



2r 
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Figure Caption: 

Figure 1: Graph of the effective potential (5.11) with (5.12) due to light pseudoscalar 
mesons. For comparison the graph of the quadratic approximation (5.13) is also shown. 
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